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Matter power spectra in dynamical Dark Energy cosmologies 
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ABSTRACT 

We used a suite of numerical cosmological simulations in order to investigate the effect 
of gas cooling and star formation on the large scale matter distribution. The simulations 
follow the formation of cosmic structures in five different Dark Energy models: the fiducial 
ACDM cosmology and four models where the Dark Energy density is allowed to have a 
non-trivial redshift evolution. Each simulation includes a variety of gas physics, ranging from 
radiative cooling to UV heating and supernova feedback (although the AGN feedback is not 
incorporated). Moreover, for each cosmology we have a control run with dark matter only, 
in order to allow a direct assessment of the effect of baryonic processes. We found that the 
power spectra of gas and stars, as well as the total matter power spectrum, are in qualitative 
agreement with the results of previous works not including the AGN effects in the framework 
of the fiducial model, although several quantitative differences exist. We used the physically 
motivated halo model in order to investigate the backreaction of gas and stars on the dark 
matter distribution, finding that it is very well reproduced by simply increasing the average 
dark matter halo concentration by 17%, irrespective of the mass. This is in agreement with 
gas cooling dragging dark matter in the very center of halos, as well as adiabatic contraction 
steepening the relative potential wells. Moving to model universes dominated by dynamical 
Dark Energy, it turns out that they introduce a specific signature on the power spectra of 
the various matter components, that is qualitatively independent of the exact cosmology 
considered. This generic shape is well captured by the halo model if we blindly consider 
the cosmology dependences of the halo mass function, bias, and concentration. However, 
the details of the dark matter power spectrum can be precisely captured only at the cost of 
a few slight modifications to the ingredients entering in the halo model. The backreaction 
of baryons onto the dark matter distribution works pretty much in the same way as in the 
reference ACDM model, in the sense that it is very well described by an increment in 
the average halo concentration. Nonetheless, this increment is less pronounced than in the 
fiducial model (only ~ 10%), in agreement with a series of other clues pointing toward the 
fact that star formation is less efficient when Dark Energy displays a dynamical evolution. 



1 INTRODUCTION 
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Many independent cosmological experiments point nowadays to- 
ward a concordance cosmological model in which gas and stars, 
that is standard baryonic matter, constitute only a few percent of 
the total energy budget of the Universe. The bulk of this budget 
is subdivided between a non-baryonic Cold Dark Matter (CDM) 
component, which decide the dynamics of cosmic structures, and a 
smooth Dark Energy (DE) component, that affects only the global 
geometry and expansion rate of the Universe, partic ularly being 



responsible for its recent accelerated expansion ph ase (Riess et al 



19981: Perlmutter et al. 1999; Tegmarketal. 2006; Komat su et al 



201 ll ; [Sullivan et al.l201 lh . The vast majority of cosmological tests 
are in agreement with a constant DE density (dubbed cosmological 



constant and indicated with A), hence the standard cosmology is 
usually labeled as ACDM. Despite its striking observational suc- 
cess, a cosmological constant is largely unappealing from the theo- 
retical point of view, for a variety of reasons that range from quan- 
tum gravity considerations to fine tuning and coinciden ce problems 
dCarroll. Press. & Turned 1 9921 ; |Peebles& Ratrall2003o . 

This fact triggered an enormous amount of work in the past 
two decades, aimed at finding alternative explanations for the 
observed geometry and accelerated expansion of the Universe 
without a cosmological constant. One such alternative consists 
in considering DE as the energy density of a scalar field (the 
quint e ssence) that slowly rolls dow n a potential well JWetterich] 
1 19951 : iHebecker & Wetterichl 1200 lh . This approach would solve 
many of the theoretical problems of a cosmological constant, how- 
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ever infinitely many such models can be devised, depending on the 
shape of the scalar field potential and on the nature of its coupling 
to gravity and matter. As it turns out, having a quintessence model 
that reproduces the recent accelerated expansion of the Universe 
within experimental errors is relatively easy, therefore signatures 
distinctive of quintessence models should be mainly searched for 
in the formation of cosmological structures. 

For the standard ACDM cosmology, the formation of dark 
matter structures is relatively well understood, thanks to both the 
linear perturbation theory and numerical cosmological simulations. 
The latter are of particular importance, in that they allow to follow 
the complete growth of structures in the deeply non-linear regime, 
for which semi-analytic methods are unsuited. The reason for this 
relative success resides in the fact that dark matter is assumed to 
respond only to gravity, which is easily modeled. Gas physics, 
on the other hand, is far more complicated, including a plethora 
of non-gravitational phenomena such as preheating, radiative and 
metal line cooling, star formation, energy feedback, etc. Many of 
these processes are active on scales too small to be resolved by cur- 
rent cosmological simulations, thus requiring the adoption of semi- 
analytic models in order to describe sub-grid physics. Despite being 
subdominant with respect to dark matter, baryonic matter can have 
a substan tial impact on the dist ribution of material at small scales 
(see, e.g.. |Puchwein et al]|2005h . and hence needs to be taken into 
account when formulating predictions for future high precision cos- 
mological tests. 

As mentioned above, given the relevance of the DE problem 
for modern cosmology, it is important to understand in detail how 
the entire process of structure formation in quintessence models 
differs from the case of a mere cosmological constant. This is- 
sue has been addressed several times in the past with the use of 
numerical simulations i Dolag et al]|2004l : iBartelman n et alj|2005l : 
iGrossi & Springel 200S ), however the focus has always been more 
on the properties of dark matter structures (abundance, distribution, 
internal features, etc.) rather than the large scale matter distribution 
in general. In addition, almost none of these works included stan- 
dard baryonic matter in the problem. Changes in the properties of 
dark matter structures are expected to affect baryonic processes, 
which depend substantially on the environment and on the inter- 
play between the cooling timescale and the dynamical timescale 
for structure collapse. At the same time, baryons are expected to 
provide some kind of backreaction on the dark and total matter dis- 
tributions, which would also depend on the background expansion 
history of the Universe. 

In this paper we thus explored in detail the co-evolution of 
baryons and dark matter in cosmologies with dynamical DE. We 
made use of the power spectra of the various matter components, as 
measured in a suite of hydro-cosmological simulations. The power 
spectrum includes a wealth of information about the matter density 
field, and it is the main subject of actual measurements such as cos- 
mological weak lensing. The simulations we used are larger than 
many, if not all, simulations used before for similar studies, includ- 
ing the mere study of baryons in the concordance ACDM cosmol- 
ogy. As a consequence, they are less affected by finite volume ef- 
fects, although they conceivably have a lower resolution than a few 
of the previous works. Large scale structure formation was evolved 
in five DE cosmologies, namely a fiducial ACDM model, two stan- 
dard quintessence models and two extended quintessence models. 
This allowed us to explore differences and regularities amongst 
models with substantially different redshift evolutions of the DE 
equation of state parameter. 

In this work we focused exclusively on the simulation snap- 



shots at z = 0. As will be seen, our set of results is already signifi- 
cantly dense in that case, thus we believe that a study of the redshift 
evolution of the matter power spectra would be better suited for a 
stand-alone effort. The rest of this work is organized as follows. 
In Section [2] we summarize the current state of knowledge about 
how the matter correlation function is affected by gas physics on 
one hand, and by a background dominated by dynamical DE on the 
other. This will help to set the stage for what follows, and to better 
motivate our work. In Section[3]we describe the five DE models (a 
fiducial ACDM and four dynamical DE) that have been explored 
in this work, and in Section[4]we give some details about the simu- 
lations employed. In Section[5]we describe the halo model, a flex- 
ible and physically motivated recipe for computing the dark mat- 
ter power spectrum that we subsequently used for interpreting our 
findings. In Section|6]we present our results, while in Section|7]we 
summarize our conclusions. 



2 PREVIOUS STUDIES 

In this Section we summarized the results of previous studies in 
the field that we are investigating, in order to better highlight the 
elements of novelty present in our analysis. 



2.1 Gas physics 

The study of the effect of baryons on the large scale matt er dis- 
tributi on in a ACDM cosmolo gy dates back at least to I White! 
(2004) and Zhan & Knox (2004). In the first paper, the author es- 
timated the changes in the total matter power spectrum due to 
cooling mechanisms that would lead the gas to concentrate in the 
very center of dark matter halos, thus dragging along dark matter 
due to the steepened potential wells. This would lead to an effec- 
tive increase in the concentration of cosmic structures, translating 
into an increase in t he matter correlation amplitude at k > 5 h 
(2004) us ed a simple model of adiab atic contraction 
dBlumenfhal et alJI 19861 : iKochanek & Whitell200lh in combination 
with the halo model (see Section|5]below) in order to replicate this 
effect, and estimated an increase in total matter power compared 
to the dark matter only case of ~ 5 - 20% at k > 10 h Mpc~', 
depen ding on the values of his model parameters. IZhan & Knoxl 
(2004) also employed the halo model, suitably modified in order 
to include the baryonic contribution to the density distribution of 
cosmic structures. In this case the authors ignored the baryons that 
cool down and condense at the center of dark matter halos, while 
focusing on the hot diffuse baryons, hence finding a decrement of 
the total matter power spectrum as compared to the dark matter 
only power spectrum. 

Numerical simulations of structure formation have been self- 
co nsistently used in order to investigate this issue for the first time 
bv lJing et al. (2006). Their simulations followed the evolution of 
2 x 512 3 gas and dark matter particles within a box of 100 hT 1 
Mpc, including radiative cooling, star formation, and supernova 
feedback. The authors found that the total matter power spectrum 
was suppressed with respect to the dark matter only case at in- 
termediate scales (k ~ I — 2 h Mpc -1 ) by a few percent, and 
then very much increased at very small scales (k > 10 h Mpc -1 ). 
This result can be se en as a comb ination of the effects of cooling 
baryons described by [jVhite (2004) and of hot baryons described by 
Zhan & Knoxl f2004). Moreover, they were able to study the back- 
reaction of baryons on the dark matter power spectrum, showing 
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that the latter is enhanced at intermediate/small scales with respect 
to the dark matter o nly case, e.g., by ~ 10% at k ~ 10 h Mpc -1 

Subsequently, Rudd, Zentner, & Kravtsov ( 2008) performed 
an analysis similar to the one of Ijing et alj J2006I) by adopting 
2 x 256 3 dark matter and gas particles in a 60 /i -1 Mpc side box, in- 
cluding gas cooling, star formation, supernova feedback, metal en- 
richment, and UV heating in their sim ulations. They fo und results 
in qualitative agreement with those of Ijing et alj J2006h . although 
several noticeable quantitative differences exist. Most remarkably, 
the raise in the clustering strength of the dark matter and total mat- 
ter at small scales, as compared to the dark matter only spectrum, is 
much more enhanced, leading to a ~ 50 - 70% increase at k ~ lOh 
Mpc -1 . The authors attribute this difference to the different imple- 
mentation of gas physics and to finite volume effects, their simula- 
tion box being rather small. 

By looking for a physical interpretation of the small scale 
enhancement in the dark matter a nd total matter power spectra, 
Rudd, Zentner, & Kravtsov! d2008h noted that dark matter halos 
tend to be ~ 10-20% more concentrated when gas cooling and star 
formation are included, in a way that is not very much dependent 
on the mass. Also, the concentration of the total matter distribution 
inside cosmic structures tends to increase more dramatically, and 
more so for smaller mass halos, where gas cooling is expected to be 
more efficient. Hence, motivated by the halo model, they corrected 
the dark matter only power spectrum by the Fourier transform of a 
representative dark matter halo density profile having concentration 
augmented by ~ 70%. They found that this correction reproduces 
fairly well the total matter power spectrum measured in their sim- 
ulations. The halo model interpretati on of the effect of baryoni c 
physics was also later resumed by IZentner. Rudd, & Hul (2008), 
that forecasted the constraints given by tomographic weak lensing 
on the mass-concentration relation of cosmic structures including 
gas cooling. 

More recently, Guillet, Tevssier. & Colombil feOlCh per- 
formed yet another similar study by using an even smaller 
simulation box (50 h~ [ Mpc on a side), and including radia- 
tive and metal line cooling, UV background, and star forma- 
tion with subsequent supernova feedback and metal enrichment. 
They found as well a substantial increase in the total matter clus- 
tering strength fo r k > 10 h Mpc -1 , indeed even larger than 
the one found by iRudd. Zentner. & Kravtsov! (2008). In order to 
mimick this, they not only increased the concentration of dark 
matter halos, but also included an additional component to the 
halo model, representing the central baryon condensation . In- 
terestingly enough, bo th Rudd, Zentner, & Kravtsov! J2008I) and 
iGuillet. Tevssier. & ColombU feOlCh did not find the suppression 
in power of the total matte r at intermediate scales that was found 
instead bv ljing et ail d2006t). 

Finally, van Daa len et alj feOllh performed a similar study 
based on numerical simulations, adopting 2 x 512 3 gas and dark 
matter particles within a box of 100 h' 1 Mpc on a side, and in- 
troducing (for the first time) the effect of energy feedback from 
Active Galactic Nuclei (AGN henceforth) to the physics imple- 
mented by earlier simulations. First of all, they found a percent 
level decrease in the clustering amplitude of total matter at in- 
termediate scales , and a subseq uent increase at small scal e s, like 
Ijing et all J2006h and unlike [R udd, Zent ner. & Kravtsov! (2008) 
and lGuillet. Tevssier. & Colombil d2010h . They also found that the 
decrement is much enhanced by AGN feedback, since more gas 
is pushed to large scales, and since star formation is suppressed, 
the increa s e at sm all scales is also suppressed. Finally, similarly to 
Ijing et al.1 (2006), they studied the backreaction of gas onto dark 



matter finding, in the absence of AGN feedback, a monotonic in- 
crease in the dark matter clustering amplitude, up to ~ 10% at 
k~ 10 A Mpc -1 . 

A study of the impact of bar yonic physics on the matter power 
spectrum was also presented in Casari ni et alj d201ll) (2 x 256 3 gas 
and dark matter particles within a 256 hT 1 Mpc side box, modeling 
radiative cooling, star formation, supernova feedback and UV back- 
ground). They did not find any reduction in the total matter corre - 
lation strength at intermediate scales, contrary to ljing et"al] j2006t) . 
Moreover, the increase at small scales of the total matter power 
spectrum was ~ 40% at k ~ 10 h Mpc -1 , intermediate between the 
results of Jing et all d2006l) and those o f Rudd, Ze ntner. & Kravtsov! 
fe008h : [Guillet. Tevssier. & Colombil d2010t) . 



2.2 Dynamical DE 

One of the most recent studies concerning how a quintessence- 
dominated background affects the large sc ale matter dist r ibutio n 
in the Universe is the already mentioned ICasarini et al. I d201ll) . 
There the authors studied the matter power s pectra in dynamical DE 
cosmologies, employing previo u s results (Linder & White! 1 20051: 
[ Francis, Lewis. & LindeJ 1 20071 ; ICasarini. Maccio, & Bonometto 
2009). Namely, the non-linear matter power spectrum in a cos- 
mological model with dynamical evolution of DE can be obtained 
from the non-linear matter power spectrum of a model with a con- 
stant w (although depending on the redshift of interest). This phe- 
nomenological recipe gives results that are accurate at better than 
~ 1%, for the dynamical DE models where it has been tested. The 
authors found that in a cosmological model where the DE equa- 
tion of state parameter evol ves according to the simple recipe of 
IChevallier & Polarskil d20oTh . the gas power spectrum is increased 
at all scales, especially at very small scales, k ~ 10 h Mpc -1 . The 
dark matter only spectrum is also increased at all scales, although 
to a lesser extent, and the total matter power spectrum is slightly 
increased at intermediate scales (~ a few h Mpc -1 ) and slightly 
decreased at small scales. Unfortunately, the authors do not go to 
a great length in discussing the physical interpretation of their re- 
sults, being mainly interested in showing that their phenomenolog- 
ical recipe works also with baryonic matter components. 

Earlier efforts performed in order to understand the impact of 
a time varying equation of state parameter for DE (we exclude the 
studies dealing only with a constant w, or where the quintessence 
field is coupled to matter) on the large scale matter distribution 
never included the effects of baryons, and date back to Klyp in et al.1 
d2003l) . They performed simula tions with varied box sizes and 
mass resolutions, considering the Ratra & Peebles] dl988l) standard 
quintessence model (see Section[3]below) with a erg matching the 
value for the fiducial ACDM cosmology, finding a substantial in- 
crease in power at s cales k> 5 h Mpc -1 , especially at high redshift. 

More recently, Ma ( 2007) performed additional «-body simu- 
lations using sever al different recipes for w(z), in cluding the simple 
parametrization of Chevallier & Polarski (2001). She finds qualita- 
tively different results from I Klypin et alj d2003l) . in that the ratio 
of the dark matter power spectrum to its counterpart in the fidu- 
cial ACDM cosmology has a depression at mildly non-linear scales 
k ~ 0.5 - 1 h Mpc -1 and a subsequent increase at small scales. The 
difference with respect to the previous work might possibly reside 
in the fact that she lets the power spectrum normalization cr 8 be 
lower in the various DE models, in order to compensate for the en- 
hanced growth factor. As a matter of fact, when the author imposes 
the linear matter power spectrum at z = to be the same as in the 
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fiducial ACDM cosmology, her results are in more agreement with 
those of iKlypin et all d2003b . 

Finally, in the aforementioned Francis, Lewis, & Linder 



2007) the authors considered again the same Chevallier & Polarski 



2001) parametrization of the DE equation of state parameter (their 
simulations have a 256 h' 1 Mpc box side length with 256 3 par- 
ticles) and, matching the cr 8 and the angular diameter distance to 
the last-scattering surface, found a sub-percent level deviation from 
the fiducial ACDM case at intermediat e scales k ~ 3 h Mpc~' and 
z = 0. These results are dissimilar from Ma ( 2007) , and more in ac- 
cordance, at least qualitatively, wifh lKlypin et alj 20031) . The quan- 
titative differences with respect to the latter can be ascribed to the 
different w{z) that have been used. 

In summary, the effect of baryonic physics on the large scale 
matter distribution has been studied in several works, employing 
both numerical simulations and semi-analytic modeling. The re- 
sults broadly agree at the qualitative level, however they display 
substantial quantitative differences. Our simulations, described in 
Section [4] have a larger box size than many of the works detailed 
above and are hence less prone to finite volume effects. For a few 
of these works we also employ a better mass resolution, and can 
thus better resolve small scale physics such as the star formation 
process. This can be useful in solving, e.g., the ambiguity about the 
slight drop at intermediate scales in the total matter power spec- 
trum. The impact of dynamical DE on structure formation is much 
less established. First of a ll, many studies used the ve ry simple 
parametrization for w(z) of IChevallier & Polarskil ( 1200 lh . which is 
not physically satisfactory since it produces an infinite DE density 
at high redshift unless the DE itself crosses the phantom barrier. 
Second, the results depend substantially on how the other cosmo- 
logical parameters are set. For instance dynamical DE models with 
the same erg as the fiducial ACDM cos mology are likely to spoil 
the a mplitude of the CMB fluctuations dFrancis. Lewis. & LindeJ 
l2007h . In the present work we used a set of physically motivated 
quinte ssence cosmologies with substantially different growth his- 
tories (De Boni et al . 2010). In all cases the erg value was adjusted 
so that the CMB fluctuations amplitude is conserved at the value 
measured by WMAP. Finally, and most importantly, in this paper 
we studied precisely the co-evolution of dark matter and baryons 
in dynamical DE cosmologies, a goal that has never been accom- 
plished at this level of detail and with such an accurate modeling of 
the gas physics. 



3 COSMOLOGICAL MODELS 

The simulations we investigated in this work span a wide range 
of flat DE cosmologies. As a reference fiducial model we adopted 
the standard ACDM model, wit h parameters in ag reement with the 
WMAP- 3 years data release iSpe rgelet all 12007). Accordingly, 
the matter density parameter is Cl m = 0.268, the dark energy den- 
sity parameter is £1 A = 1 - fi nl _o, and the baryon density param- 
eter is f2t,,o = 0.044. The Hubble constant reads Hq = /il00 km 
s -1 Mpc -1 , with h = 0.704. Finally, the primordial power spec- 
trum normalization is set by cr 8 = 0.776, while its slope is fixed 
to ns = 0.947. In what follows, quantities related to this fiducial 
model will always be labeled by the suffix F. The main difference 
between the cosmological interpretation of the WMAP-3 year data 
release and those of subsequent releases is in the value of the power 
spectrum normalization, which tends to be relatively larger for the 
latter as compared to the former. However we note that the value of 
erg adopted here is compatible at 68% Confidence Level with the 



value extracted from the latest WMAP-7 release dKomatsu et al.1 
1201 lb . 

We then considered four different models involving a dynam- 
ical evolution of the DE component. For all these models the cos- 
mological parameters are the same as in the fiducial scenario, ex- 
cept for <r 8 that was rescaled in a way such that the linear mat- 
ter power spectrum is always CMB-normalized. Specifically, if 
D+(z) is the linear growth factor of a specific cosmology, and 
g(z) = (1 +z)D + (z), then 



0"8 = 0"8f' 



gFfecMB) 
g(ZCMB) 



(1) 



where zcmb = 1089. It is important to note that, since the growth 
factor in quintess ence cosmologies is always larger than its fiducial 
counterpart (see lDe Boni et al.ll2010b . the power spectrum normal- 
ization will always be such that tr 8 < cr 8)F . 

We now proceed to present in more detail the four dynami- 
cal DE cosmologies that we investigated. The first two are standard 
quintessence models, in which the DE is just the energy density as- 
sociated wi th a minimally coupl e d scalar field rol ling down some 
potential dRatra & Peebleslll988l : IWetterich!ll988b . The other two 
instead arise from a non-minimal coupling of the scalar field with 
gravity, giving effectively origin to the so-called scalar-tensor grav- 
ity theories. In what follows we briefly describe both pairs of mod- 
els separately. 

3.1 Standard quintessence 

In standard quintessence models the evolution of the scalar field 
<p whose energy density represents DE is governed by the Klein- 
Gordon equation (we work in natural units, c = h = 1 unless other- 
wise noted) 



if + 3H<p + 



dV 
dtf 



0. 



(2) 



which is a direct consequence of the Friedmann equations. In Eq. 
(|2j a dot denotes derivative with respect to cosmic time, while V(ip) 
is the potential of the scalar field. 

As can be seen, the evolution of the scalar field if, is deter- 
mined by the expansion rate of the Universe, which in turn depends 
on its whole matter-energy content, and by the effective potential of 
the scalar field. The two standard quintessence models considered 
in this work arise from two different choices for the dependence 
of V on if. In the first case a power-law dependence is adopted, 
according to which 



V( V ) = A* + "— a 
tf a 



(3) 



dRatra & Peebles 1988b , where both X and a ^ are free parame- 
ters. This model will be labeled as RP in the remainder of this work. 
The second model arises from a generaliz ation of the RP poten tial, 
suggested by supergravity considerations dBrax & Martiij2o"o(il) . In 
this case the potential reads 



vbp) = r 



'i;exp(4^ 2 ), 



(4) 



and will be labeled as SUGRA henceforth. 



3.2 Scalar-tensor theories 

For the second class of models the scalar field tp is coupled non- 
minimally with gravity, giving rise to an effective scalar-tensor 
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gravity theory. We call these models extended quintessence models. 
In order to be more specific, a new function f(tp,K) is to be intro- 
duced in the gravity action, which describes the coupling between 
the quintessence scalar field and the Ricci scalar R. In the extended 
quintessence cosmologies considered in this paper, we assume that 
the coupling function can be factorized as 

f(.<p,R) = R[l+&nG,t(<p 2 -<pl)]. (5) 

In the previous equation, G, is a 'bare' gravitational constant, that 
is in general different from Newton constant G and it is set in order 
to match local constraints on General Relativity. 

The strength of the coupling is defined by the parameter 
Particularly, in this paper we have considered two extended 
quintessence models, one with £ < (specifically £ = -0.072, EQn 
henceforth), and one with £ > (£ = +0.085, EQp). The potential 
of the scalar field, although not the only source term of the resulting 
Klein-Gordon equation, needs to be specified in this case as well. 
Here, we adopted the RP power-law potential described above. It 
is worth noticing that, while the 'bare' gravitational constant G, is 
a constant of the theory, the effective gravitational constant that is 
experienced by particles in the simulation is a function of cosmic 
time. Specifically, it converges to G, at z = and approaches a 
constant value larger (smaller) than G, at high redshift for positive 
(negative) £ . 



having a mass of 4.4 x 10 9 M Q h- 1 . In the following, we will often 
refer to the dark matter power spectrum in a control run as simply 
dark matter only power spectrum. 

The power spectra of each matter component in each simu- 
lation were computed by projecting the matter distributions on a 
regular grid constituted by n = 960 grid points on a side, and then 
using discrete Fourier transforms. The matter density at every grid 
point wa s evaluated by the Triangular Shape Cloud (TSC) method 
(refer to iHocknev & Eastwood Il988h . which was then compen- 
sated for interpolation effects (see also |jing||2005h . The Nyquist 
frequency corresponds therefore to £ N = (n - 1) n/L ~ lOh Mpc~', 
meaning that the matter distribution is not sufficiently sampled 
(and hence was not considered in our analysis) on scales k > fc N . 
The simulations we used were not re-run adopting mass resolu- 
tions different from the nominal one, and hence we have no way to 
straightforwardly asses the robustness of our power spectra analysis 
with respect t o partic le number and/or box size. Nonetheless, e.g., 
iBorgani et al.l d2006l) shown that changes in the mass resolution 
have a negligible impact on the simulated distribution of baryons 
at cluster scales, except for a change in the star formation rate at 
z > 3, thus hinting that this should not be an issue for the present 
work. 



4 NUMERICAL SIMULATIONS 

The n umerical simulations are fully described in |Pe Boni et al.l 
(2010). Here we just recall the main features, referring the 
interested reader to that work for further details. The evo- 
lution of the large-scale matter distribution in the various 
DE cosmologies was simulated by using the GADG ET- 3 code 
dSpringel. Yoshida. & Whitel200lHSpringel et al]200lh . This code 
allows to simulate the evolution of both the dark matter and the gas 
comp onent, making use of the en tropy conserving formulation of 
SPH dSpringel & Hernquistl2002h . The code was suitably modified 
in order to allow for an ar bitrary time evolut ion of the DE equation 
of state (see for instance [Dolag et"alj 2004), an d also for an arbi- 
trary time evolution of the effective gravitational constant, needed 
for the extended quintessence models. 

The baryonic physics implemented in the simulations includes 
radiative coolin g, heating by an un iform redshift-dependent UV 
background (Haardt & Madau 1996), and star formation and feed- 
back. The latter requires sub-grid modeling of the multi-phase Inter 
Stellar Medium (ISM), with the cold phase providing the reservoir 
for the star formation, and the hot phase being fueled by supernova 
heating. The latter provides also evaporation of some cold clouds, 
thereby leading to a self-regulation of the star formation process. 

Each simulation has a comoving box side measuring L = 
300 lr [ Mpc and follows the evolution of 768 3 dark matter particles 
and an equal number of gas particles. The mass of each dark matter 
particle is 3.7 x 10 9 A f o fe~ 1 , while that of each gas particle equals 
7.3 x 10 s M /T>. As in lDolag et all d2004h . the initial redshift is not 
the same for each simulation, rather for each DE cosmology it is 
rescaled in order to conserve the linear growth factor at the initial 
redshift in the fiducial ACDM simulation. Accordingly, all simu- 
lations start from the same random phases, but with amplitude of 
the initial fluctuations rescaled to satisfy the CMB constraints. For 
every cosmological model, a dark matter only control simulation 
was ran, providing the reference for the effect of gas physics. In the 
control runs the number of dark matter particles is 768 3 , each one 



5 THE HALO MODEL 

In the course of our analysis we wanted to have a semi-analytic pre- 
scription for computing the dark matter power spectrum, in order 
to physically interpret our re sults. While fits to numerically sim- 
ulated matter power spectra jPeacock & Do dds 1994: [Smith et al.l 
120031) are widely used, they have two drawbacks for our pur- 
pos es. First, they do n ot really have a physical foundation (although 
the lSmithetalj|2003l prescription uses some elements of the halo 
model, it is ultimately only a fit), thus they are not useful when try- 
ing to interpret our results in a physical framework. Second, these 
fits have been matched against simulations containing dark mat- 
ter only and in a ACDM background. Since they lack a physical 
motivation, there is no telling what i s their ran ge of validity for al- 
ternative cosmologies (see however Ma 2007), or even for scales 
outside the range probed by those simulations. As a consequence, 

we adopted the physically m otivated hal o model . 

The halo mod el dMa & Frvl |2000| : ISeliakl l2OO0l: 
ICoorav & Shetbl |2002|) is a physically motivated framework 
for describing the large scale matter distribution. It is based on 
the assumption that each matter particle in the Universe is locked 
inside halos of some mass. If this is true, then the correlation 
function of matter particles is given by the sum of two contribu- 
tions. The first one is the contribution of particle pairs belonging 
to the same halo, and dominates on small scales. The second one 
is the contribution of particle pairs residing in separated halos, and 
dominates on large scales. It is natural to expect that, while the 
2-halo term depends on the bias of dark matter halos, the 1-halo 
term is more dominated by the internal structure of individual 
halos. For details on the imple mentation of the hal o mod el we 
refer the interested reader to iFedeli & Moscardirul d2010l) and 
references therein. Here we just want to stress the particular 
importance of the relation between mass and average concentration 
of cosmic structures. In the present work we adopted an empirical 
relation that reproduces well the power spectrum fits to n umerical 
simulations dPeacock & Doddsl 1 1994 ISmith et al.l 120031) for the 
ACDM cosmology, that is 
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MM 



(6) 



where Co = 10 and f} = -0.15. In the previous equation M,(z) is 
the typical collapsing mass at a given redshift, and it is the only 
cosmology-dependent quantity that enters the concentration-mass 
relation. This means that when we applied the halo model to cos- 
mologies with dynamical DE, the concentration-mass relation for 
dark matter halos was modified only through Af,(z). We also ob- 
serve that the concentration of halos with a given mass decreases 
faster with red shift according toEq. j6) than the oc (1+^) -1 behavior 
found by, e.g.. lBullock et al.l J200lh or lEke. Navarro. & Steinmet3 
d200lh - This kind of issue has been already noted and studied in 
detail bv lHuffenberger & Seliakl J2003h . 

Before proceeding we would like to stress a subtlety of our ap- 
proach. In the present work we tried to fit ratios between simulated 
power spectra with ratios between outcomes of the halo model. Al- 
though the halo model by itself gives a good description of the true 
matter power spectrum, it cannot be expected to be precise at the 
percent level, as are the effects that we are considering here. This 
is due to the many obvious approximations and simplifications that 
are introduced in the model. Besides, the prescriptions for the mass 
function of dark matter halos and their linear bias (that are two 
ingredients of the halo model) can be considered accurate only at 
the ~ 10% level at best. However, by considering only ratios we are 
presumably canceling out these discrepancies, thus highlighting the 
true physical effects. 



6 RESULTS 

We now proceed to present our results. We treat at first the effect 
of baryonic physics in the fiducial ACDM cosmology, in order to 
compare our results with those of previous works in the field. Then 
we move on to consider the effect of a dynamical DE evolution 
on each matter component of the simulations, including both dark 
matter and baryons. Subsequently we explore the impact of bary- 
onic physics on the dark matter distribution in each quintessence 
cosmology, and we conclude by exploring the behavior of the total 
matter power spectrum in these models. 

In what follows, we shall always refer to the dimensionless 
power A 2 (k,z), defined as 



2tt 2 



(7) 



The dimensionless power is basically the contribution to the mass 
variance given by logarithmic intervals in wavenumber. As men- 
tioned, in many circumstances we have shown ratios of powers, 
thus rendering the distinction between the dimensionless power and 
the power spectrum irrelevant. 

6.1 Impact of gas physics in the ACDM cosmology 

We first turn attention to the power spectrum of the various matter 
components in the ACDM simulation. For definiteness, in FigureQ] 
we show the z = dark matter only power for the fiducial ACDM 
cosmology control run, as well as dimensionless powers for the 
different matter components in the full simulation. 

The most evident fact about this Figure is the large discrep- 
ancy between the power spectrum of the stars and that of the other 
matter components. This difference is expected, since stars tend to 
form inside dark matter halos, that are biased tracers of the un- 
derlying smooth density field. Therefore stars are more clustered 
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Figure 1. The dimensionless powers of dark matter, gas, and stars in the 
fiducial ACDM cosmology. The black solid line refers to the dark matter- 
only simulation, while the other line types refer to the various matter com- 
ponents in the full simulation. We did not plot the total dimensionless power 
since on these scales it basically overlaps with the dark matter one. All 
curves refer to z = 0. 



together compared with dark matter or gas particles. A similar be- 
havior for the correlati on function of sta rs in numerical simulations 
has been found, e.g., in Jing etal. (2006) (see Section[2]above). 

Next, we consider the behavior of gas. As can be seen by look- 
ing at Figure [T] the gas traces quite well the dark matter distribu- 
tion at large scales, with only very little (and smaller than unity) 
bias. However, as the scales become more and more non-linear, the 
clustering strength of gas particles becomes smaller and smaller as 
compared to that of dark matter. This behavior is due to a combi- 
nation of star formation and hot gas pressure. Specifically, part of 
the gas tends to cool down and transform into stars in the cores of 
dark matter halos, implying a depletion at relatively small scales. At 
the same time, shock heating and gas pressure prevent the hot gas 
component to have small scale clustering. Agai n, similar qualita - 
tive conclusions have been reached prev iously by Jing et al. d2006h . 
iRudd, Zentner, & Kravtsovl d2008h . and lvan Daalen et al.l l l201lh . 

Figure Q] also allows us to understand what happens to the 
dark matter component upon inclusion of the gas physics, that is 
the backreaction of gas upon dark matter, which is a central fo- 
cus of this work. As can be seen, the dark matter indeed gains a 
little bit of clustering strength at small scales. This can be inter- 
preted in the following way: while the gas cools and concentrates 
in the very center of dark matter halos, the dark matter particles 
are also pulled closer to the halo centers. Moreover, adiabatic con- 
traction due to gas radiative losses also contributes to make dark 
matter halos more concentrated. Overall, dark matter particles in- 
side individual structures tend hence to be m ore strongly clu stered. 
The same qualita t ive ef fect has been found bv ljing et al.ld2006l) and 
Ivan Daalen et alj d201 lh . We can further investigate this important 
issue by plotting the ratio of the various spectra to the dark matter 
only spectrum. This is done in Figure|2] 

In this Figure we did not show the stellar power spectrum, 
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Figure 2. The ratio between the power spectra of the various matter com- 
ponents in the full simulation (as labeled in the plot) for the fiducial ACDM 
cosmology and the matter power spectrum for the dark matter only con- 
trol ran. The yellow circles represent the halo model predicted dark matter 
power spectrum obtained by increasing the average structure concentration 
according to the recipe detailed in the text, normalized by the original halo 
model spectrum. All quantities refer to z = 0. 



since it would be way out of scale. On the other hand, we added 
the total matter power spectrum, that in Figure Q] would be too 
similar to the dark matter one to be visible. In this plot we can 
better appreciate the different behavior of the gas component as 
compared to the dark matter one. Particularly, the clustering am- 
plitude of the former is suppressed by up to a factor of ~ 2.5 in 
comparison to the latter at scales k < 10 h Mpc -1 . In several previ- 
ous works dRudd Zentner. & Kravtsov||2008l ; ICasarini etal]|201ll ; 
Ivan Daalen et al.ll201 ll) it can be seen that the gas power spectrum 
has the tendency to catch up on the dark matter-only spectrum at 
scales k ~ 10 h Mpc~' or smaller. We do not see trace of this be- 
havior in our plot, although it is possible to happen at scales below 
the smallest ones probed by our runs. We attribute this discrepancy 
to the di fferences in the implementation of gas physics. Finally, as 
noted bv lRudd. Zentner. & Kravtsovl (I2008). the large scale smaller 
than unity bias of the gas component with respect to the dark matter 
component is likely due to the fact that low mass dark matter halos 
(that are less biased with respect to the large scale matter distri- 
bution) contain a larger gas fraction as compared to the high mass 
(and more clustered) structures. 

Let us now focus on the modification to the dark matter power 
spectrum due to the inclusion of gas physics, that is the backre- 
action. As can be seen from Figure [2] the dark matter power is 
correspondingly enhanced by up to a factor of ~ 15 - 20% at 
the smallest scales probed by the simulatio ns. This is somewhat 
larger than what is found bv ljing et al.1 d2006h and lvan Daalen et all 
d201ll) (~ 10%), and at the same time smaller than the results of 
ICasarini et al. I Hoi 1) (~ 40 - 50%). Again, these differences can 
be explained in t he different details of gas physics implementation 
that, as shown by van Daalen et al] d201ll) can have a large impact. 
In Figure [2] we also plot the dark matter power spectrum evaluated 



through the halo model by assuming a higher average concentra- 
tion of dark matter halos, normalized by the regular halo model 
output. More specifically, we modified the concentration of a halo 
with given mass and redshift according to cp(M,z) = fcp$(M,z), 
where / = 1.17. As can be seen, this very simple recipe accounts 
very well for the effect observed in the simulations, and the devia- 
tions between the halo model prediction and the simulation results 
are < 1% at all scales probed by the latter (< 10 h Mpc -1 ). 

Realistically, one might expect the effect of baryonic physics 
to act differently on dark matter halos of different masses, with less 
massive halos being more affected as compared to more massive 
structures. Nevertheless, simply changing the concentration of all 
halos of the same amount irrespective of the mass is much sim- 
pler and, as can be seen, works remarkably well. The modifica- 
tion to the halo model described above matches also fairly well the 
total matter power spectrum, represented by the cyan line in Fig- 
ure [2] with maxim um deviations being of the order o f 3 - 4%. As 
a matter of fact in lRudd. Zentner. & Kravtsovl J2008t) . the authors 
used halo model considerations and increased the concentration of 
st ructures in order to reproduce the to tal matter power spectrum. 
In lGuillet. Tevssier. & Colombil d2010h the authors did the same by 
also including a model for the distribution of stars in the very center 
of dark matter halos. Both works managed to reproduce the ratio of 
the total matter power spectrum to the dark matter only power spec- 
trum with an accuracy of ~ 5% or worse, and deemed this result as 
acceptable. 

It would be interesting if the halo model could be straightfor- 
wardly modified in order to reproduce the total matter power spec- 
trum at the sub-percent level of accuracy since, for instance, cosmic 
shear surveys measure the total power, rather than the dark matter 
power. However the behavior of the former is more complicated 
than that of the latter. First of all, there is a slight suppression of 
power, at the percent level, at interm ediate scales. As mentioned 
in Sec tion [2] this was not found by iRudd. Zentner. & Kravtsovl 
(2008) and Guill et. Tevssier. & Colombil ||201C|). but was found~b y 



Ijing et al. I J2006h and, more recently, by van Daal en et al. I( l201ll) . 
Even if this is ignored however, the total matter power spectrum 
tends anyway to be shallower than the dark matter one at interme- 
diate scales and steeper at small scales. This is due to the fact that 
the gas contribution pushes it down at intermediate scales and the 
stellar contribution pulls it up at small scales. As a consequence, 
the behavior of the total matter power spectrum cannot be easily re- 
produced with a simple modification of the concentration of struc- 
tures, at least if we want sub-percent level precision. From the halo 
model point of view, this can be interpreted as the consequence 
of a twofold effect: (f) the shape of individual structures (contain- 
ing dark matter, gas, and stars) is not as well described by a NFW 
profile as are dark matter halos alone; (if) the collisional nature of 
gas makes sure that the assumption according to which all mat- 
ter is locked into structures is not valid. In other words, the halo 
model should include an additional component describing the frac- 
tion of gas that is still smoothly distri buted, similarly t o what hap- 
pens for Warm Dark Matter ( ISmith & Markovidl201lh scenarios. 
An interesting direction for future investigation would be to de- 
velop a generalized halo model that would include the distribution 
of dark matter, hot and cold gas, and stars. The parameters of such 
a model could definitely be adjusted to reproduce the total matter 
power spectrum at the sub-percent level, however it is well beyond 
the scope of the present paper. 
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Figure 3. The ratio between the power spectra of the various matter components in the four quintessence simulations adopted in this work and the corresponding 
quantities in the fiducial ACDM run. The black solid line refers to the power spectrum of the dark matter only simulation, while the other line types refer to 
the various matter components of the full simulations, as labeled in the plot. All quantities are evaluated at z = 0. 



6.2 Matter power spectra in quintessence cosmologies 

In Figure[3]we show the power spectra of the different matter com- 
ponents in the four dynamical DE simulations considered in this 
work, normalized by the corresponding quantities in the fiducial 
ACDM run. Due to the fact that the normalizations of the z = 
power spectra in the quintessence models are always smaller than 
the corresponding quantity in the fiducial model, the clustering 
strength of all matter components at very large scales is always 
reduced. The first interesting thing to note is that all matter com- 
ponents have the same generic shape, being flat at large scales, as 
expected, with a depression at k ~ 1 h Mpc -1 and then growing 
again at smaller scales. This aspect is reminescent of the dark mat- 
ter pow er spectrum shape in dynamical DE models explored bv lMal 
( 2007), with the difference that there the depression was located at 
slightly larger scales. As a matter of fact, even here we note that 
the dip is at slightly smaller scales in the models RP and EQp as 
compared to the models SUGRA and EQn, suggesting that while 
the presence of the depression is a generic feature, its precise lo- 
cation depends on the exact background expansion history of the 



Universe and on the amplitude of primordial density fluctuations. 
This is supported by the fact that, as we will discuss shortly with the 
help of the halo model, the scale of the dip depends on the scale at 
which non-linear effects become important. It should also be noted 
that when the value of tr 8 is adjusted to be the same for all models 
the depression is expected to dis app ear, as happen in the r esults of 
iFrancis. Lewis. & Linderj d2007h and lCasarini et alj d201ll) . 

Let us next examine the behavior of the gas component. The 
relevant power spectrum ratio is equal to the ratio of the dark matter 
power spectra (both in the dark matter only and in the full simula- 
tions) at large scales. However, at intermediate/small scales it tends 
to be larger. This means that, while the gas clustering amplitude in 
quintessence cosmologies remains smaller than its ACDM coun- 
terpart, (0 less gas is depleted and locked into stars in cosmologies 
with a dynamical evolution of DE, and (if) shock heating and hot 
gas pressure are slightly less efficient in preventing gas clustering. 
Which one of these t wo effects is the dom inant one is hard to say, 
however we note that lDe Boni et alj d20 ldh already pointed out that 
star formation tends to be suppressed in dynamical DE cosmolo- 
gies as compared to the fiducial ACDM model. An increase in the 
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gas clustering strength at s mall scales for dynam ical DE models has 
furthermore been found bv lCasarini etal] J201 lh . where the authors 
also note that star formation is reduced. Curiously, this means that 
the reduction in the star formation efficiency is not (or not com- 
pletely) to be attributed to the smaller erg values adopted here. All 
in all this fact makes sense, because a dynamical evolution of dark 
energy generically implies a larger Hubble drag, and thus an ear- 
lier suppression of structure formation. It is plausible that this also 
translates into a suppression of the star formation process. 

The above mentioned conclusion is also supported by observ- 
ing the stellar power spectrum. Specifically, the ratio of the stellar 
power spectra is always larger than the ratio of the power spec- 
tra for the other matter components, in particular dark matter. This 
means that the clustering strength of the stars is less suppressed in 
dynamical DE simulations with respect to the clustering strength of 
dark matter. In other words, in the latter simulations the stars tend 
to form into more biased halos as compared to the fiducial run. This 
is consistent with gas cooling and star formation being more diffi- 
cult, and hence effective only in more massive halos as compared 
to the fiducial run. 

Finally, it is interesting to note that at small scales the ratio 
between the dark matter (and total matter) power spectra in the full 
simulations tends to be smaller than the ratio of the dark matter 
spectra in the dark matter only runs. This hints to the fact that the 
dark matter power spectrum is less affected by the baryonic physics 
with respect to the fiducial ACDM run, again confirming the afore- 
mentioned conclusions. Summarizing this part, we can conclude 
that evolving structures in a quintessence cosmology have an over- 
all qualitative effect on the matter distribution that does not depend 
much on the specific equation of state parameter w(z) adopted, al- 
though the details do, especially at small scales. It is also interesting 
to note that the extended quintessence models, despite having an 
effective modification of the gravitational constant, do not behave 
much differently from the others, standard quintessence models. 

In order to better understand the depression in the power spec- 
trum ratios of all matter components that is visible at k ~ 1 h Mpc~' 
for all quintessence models, we tried to model the ratios between 
the power spectra in the dark matter only simulations with dynam- 
ical DE and the corresponding quantities for the fiducial run by 
using the halo model. Preliminarily, we should recall that the halo 
model has some intrinsic dependence on cosmology, through the 
linear matter power spectrum, the halo mass function and the cor- 
responding linear bias. Moreover, the relation c(M,z) also has a 
cosmology dependence due to M,(z). All these dependences have 
been suitably taken into account, by adopting the correct cr s and 
growth factor for computing the mass variance and by solving the 
spherical collapse equations w ith quintessence contribu tions. How- 
ever, the gener ic shape of thelSheth & Tormer] J2002h mass func- 
tion and of the lSheth. Mo. & Tormenl ( 1200 ID linear halo bias were 
left intact. 

In Figure|4]we report the same as in Figure[3]but only for the 
SUGRA quintessence model and for the dark matter power spec- 
trum in the dark matter only control run. In addition, we overplot 
the prediction given by the halo model obtained by blindly cor- 
recting for cosmology only as described above. We also report the 
two contributions to the halo model power spectrum, the 1-halo 
term and the 2-halo term, again in units of the dark matter power 
spectrum in the reference ACDM cosmology. There are a few in- 
teresting features to be noted about this plot. First, the overall shape 
of the power ratio is recovered, and in particular the depression that 
is visible at k ~ 1 h Mpc~'. We can also understand the origin of 
this depression. As one might expect, the 2-halo term is a mono- 
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Figure 4. The ratio between the dark matter power spectrum in the SUGRA 
quintessence model to the corresponding quantity in the ACDM control run. 
The black solid line refers to the dark matter only simulation, as labeled. 
The yellow circles represent the halo model prediction, broken in the 1 -halo 
(cyan triangles) and 2-halo (cyan squares) contributions (see the text for 
more details). 



tonically decreasing function of the wavenumber, while the 1-halo 
term is a monotonically increasing function. Therefore, it is natural 
that the sum of the two displays a minimum, although this mini- 
mum does not correspond to the point in which the two terms have 
equal contributions because they also have a different steepness. 

Next, we note that the halo model representation is not entirely 
accurate from the quantitative point of view, for three reasons: (£) 
at very large scales the power spectrum ratio evaluated by the halo 
model is slightly larger than its simulated counterpart, by about 
~ 2%; (h) at intermediate scales this overestimation worsens, or 
in other words, the dip in the ratio is not accurately reproduced; 
(iii) at small scales the halo model prediction has a different slope 
from the simulations, so that even if the former was normalized 
down, we would still miss agreement. The same conclusions apply 
to the other three quintessence cosmologies. We now explore each 
of these issues in order, explaining how we tackled each of them 
in order to obtain a more accurate description of the dark matter 
power spectrum in dynamical DE cosmologies. 

First of all, we believe the disagreement at large scales not to 
be of physical, rather of numerical origin. At those large scales the 
ratio of the power spectra should equal the ratio of their primordial 
normalizations, that is the ratio of the respective cr$ values squared. 
While this is true for the halo model, the simulated spectra are ~ 2% 
off. The only reason for this is that the growth experienced by struc- 
tures in the simulation is different from the theoretical growth fac- 
tor. As a matter of fact, we verified that this is the case for all the 
quintessence models at the ~ 1% level, which is perfectly accept- 
able in terms of numerical precision, but that produces the observed 
offset. In order to compensate for this, we slightly shifted the theo- 
retical tr 8 values adopted, in order to get a match at large scales. 

Second, the discrepancy at intermediate scales is somewhat 
expected, since the halo model is well known not to be percent 
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Figure 5. The ratio of the dark matter power spectra in the four dynamical dark energy cosmologies studied in this work to their counterpart in the fiducial 
cosmology. The black line refer to the dark matter only run. The yellow circles represent the halo model predictions relative to the black lines, where the 
three corrections detailed in the text have been taken into account. Particularly, the changes necessary to the mass function (n/np) in order to get agreement at 
intermediate scales and to the average halo concentration (c/cp) to get agreement at small scales are labeled in each panel of the plot. 



level accurate at the transition between the 1-halo and the 2-halo 
term. It is very well possible that this accuracy actually depends on 
cosmology, and hence we see a larger discrepancy at those scales. 
There exist prescriptions that allow to improve the agreement be- 
tween the halo model and numerical simulations at intermediate 
scales, for instance the so-called halo exclusion. It means that one 
should exclude from the computation of the 2-halo term dark mat- 
ter particle pairs that are in different but overlapping halos. In other 
words, one should remove the correlations which arise on scales 
smaller than the sum of the virial radii of the two halos. At first 
approximation the halo exclusion manifests in lowering the up- 
per integration bo undary in the 2-halo term JTakada & Jain 2003; 
iTinker et alj |2005|). although more accurate prescriptions exist 
( Smit h. Scoccimarro. & Shethl2007lJSmith. Desiacques. & M arian 
l201lh . However, we verified that the halo exclusion does not lead 
to significant modifications in this case, thus implying that the halo 
exclusion itself works in quintessence cosmologies as it does in the 
fiducial ACDM model. On the other hand, we noted that the behav- 
ior of the power spectrum ratio at these intermediate scales is quite 



sensitive to the mass function. By decreasing the mass function for 
dynam ical dark energy cosmologies by only ~ 3% with respect 
to their lSheth & Tormenl d2002h predictions brings good agreement 
with the simulated curves. We therefore adopted this minor correc- 
tion, reminding that e ven for the fiducial ACD M cosmology the 
agreement between the lSheth & Tormenl j2002h mass function and 
numerical simulations is only at the ~ 10% level, and that the level 
of agreement might very well depend on cosmology. This however 
has also an interpr etation different f r om the numerical one: the nor- 
malization of the ISheth & Torrner] ( 120021) mass function ensures 
that all the mass in the Universe is locked in to halos (the basic 
assumption of the halo model). This assumption is of course not 
strictly true, as some dark matter will certainly be diffused outside 
halos. Lowering this mass function normalization by the amount 
discussed means that ~ 3% less mass in quintessence cosmologies 
is locked into halos, an interpretation that agrees with structure for- 
mation being less efficient in those models. 

Finally, the disagreement at the very small scales can be 
healed by adjusting the dark matter halo concentration, exactly as 
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Figure 6. The ratio between the power spectra of the various matter components in the full simulations (as labeled) and the dark matter power spectra in the 
dark matter only control runs. Each panel refer to a different quintessence cosmology that has been studied in this work. The yellow circles represent the halo 
model predictions assuming that the average dark matter halo concentration is increased with respect to the dark matter only case by the quantity labeled in 
each panel as c/cq. 



we did for the effect of baryonic physics. In Figure|5]we show the 
halo model predictions resulting from the three corrections outlined 
above, for all dynamical DE cosmologies considered in this work. 
In each panel we also report the corrections to the mass function 
and to the average concentration of dark matter halos that produce 
the smallest deviations. With these corrections, the power spectrum 
ratios given by the halo model reproduce their simulated counter- 
parts with a relative accuracy of < 1 - 2%. 

As can be seen by the numbers in the various panels, in order 
to bring the halo model predictions in agreement with the simula- 
tions at intermediate scales, the mass function in the quintessence 
models needs to be lowered by 2 - 4% at most, a value that is 
well within the uncertainties of the mass function even for the fidu- 
cial ACDM case. After that, in order to match the behavior of the 
power spectrum ratio at small scales, the average concentration of 
dark matter halos needs to be increased by ~ 5 - 10%, with the 
exception of the model EQp which needs no such correction. In 
any case, it is likely that the increase in concentration is not en- 
tirely physical, rather it is at least in part a response to the mass 



function decrement, that slightly suppresses the 1-halo term (not 
the 2-halo term, that is normalized at large scales). Since both ef- 
fects are needed for the overall agreement, it is hard to say what the 
actual change in concentration would physically be. 

6.3 Gas backreaction in quintessence cosmologies 

In this Subsection we address another tricky point, that is under- 
standing the effect of gas physics on the dark matter power spec- 
trum when the latter is evolved in a background cosmology dom- 
inated by dynamical DE. To that purpose, we show in Figure [6] 
the power spectra of the various matter components divided by the 
dark matter power spectrum in the dark matter only simulation for 
all the four quintessence cosmologies considered in this work. 

The behavior of all matter components looks very similar to 
that for the reference ACDM cosmology shown in Figure [2] Upon 
careful investigation however, it can be noticed that at very small 
scales the gas tends to be slightly less depleted as compared to the 
fiducial case, in agreement with the conclusions that we reached 
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before. The most interesting part of this Figure is the comparison 
of the dark matter power spectrum with the predictions of the halo 
model. In order to describe the matter power spectrum in dynamical 
DE cosmologies we adopted the modifications to the mass function 
and average halo concentrations that were found necessary in the 
previous Subsection. Then, similarly to what we have done for the 
ACDM case, we increased the concentrations in order to mimick 
the effects of baryonic physics. Since gas cooling and star forma- 
tion have been proven to be less effective when dark energy is dy- 
namically evolving, we expected to find that an increase in concen- 
tration lower than the 17% found for the fiducial cosmology would 
do the job. 

This is indeed our conclusion. With the exception of the RP 
model, that requires an increase in concentration comparable to 
the ACDM cosmology, in all other cases and excellent agreement 
can be reached with an increase of only ~ 10%, as labeled in the 
various panels of Figure [6] In all cases, the relative agreement be- 
tween the simulated dark matter power ratio and the one predicted 
by the halo model is always < 1%. In summary, the halo model 
with the modifications that we introduced given by the evolution 
in a quintessence dominated background and the effect of baryonic 
physics described above reproduces very well gas backreaction on 
the power spectrum of dark matter, even in cosmologies different 
from the fiducial one. 

The way in which the total matter power spectrum differs from 
the dark matter only one is rather similar to the fiducial ACDM cos- 
mology. Namely, also for quintessence models we can see a slight 
suppression of power at mildly non-linear scales, while the behav- 
ior at smaller scales cannot be represented by a simple change in 
the concentration. Nevertheless, if precision at the level of ~ 3-4% 
is enough, the enhanced concentration values that we quote in Fig- 
ure [6] are a good choice also for the total matter power spectra. In 
the next subsection we explore with more detail the issue of the de- 
pression in the total matter power spectrum at intermediate scales 
in quintessence cosmologies. 

6.4 The total matter power spectrum 

In this last Subsection we want to focus more on the total mat- 
ter power spectrum. As we mentioned before, it is not possible to 
obtain a precise representation of the total matter power by using 
the halo model in a way as straightforward as for the dark matter 
power. Nevertheless, understanding the behavior of the total power 
spectrum is important for weak lensing surveys. The gravitational 
deflection of light measures the large scale distribution of all mat- 
ter, irrespective of its nature. Given the large efforts that are be- 
ing dedicated to future all-sky weak lensing surveys such as Euclid 
( ILaureiisll2009n and WFIRST, it is worth studying in more detail 
this important issue. The total matter power spectra presented in 
this Subsection, complemented with their redshift evolutions, can 
be used in order to study the effect of baryonic physics in these four 
quintessence cosmologies on weak lensing observables. 

Let us begin with the fiducial ACDM cosmology. In Figure UJ 
we show, for this cosmological model, the ratio between the total 
matter power spectrum and the dark matter power spectrum in the 
dark matter only simulation. Basically we are showing the same 
cyan line that is visible in Figure [2] however this time we report 
it in percentage, in order to better highlight small differences. As 
noted above, there is a ~ 1% depression at scales k ~ 1 h Mpc~'. 
This is in good agree ment, both qualitatively and quantitatively, 
with Ijing et al.l <|2006). They found that this depression was sta- 
ble also upon exclusion of the radiative cooling and star formation 
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Figure 7. The total matter power spectrum divided by the dark matter power 
spectrum in the dark matter only simulation (cyan dot-dashed line). For 
reference, we also show the dark matter only power spectrum (black solid 
line). This panel refers to z = only. Note that the numbers on the ordinate 
axis express a percentage. 



in the gas simulation, namely by treating the gas as adiabatic. This 
mean that this feature is not related to the cold gas phase, that con- 
densates and produces stars, rather it should be related with the hot 
phase. As a matter of fact, the hot phase extends up to intermediate 
scales, where the effect of adiabatic contraction and of the stellar 
concentrations are not relevant. There, shock heating and the ther- 
mal pressure prevent the gas, and therefore the total matter, from 
strongly clustering. 

The very small amount of the depression in total cluster- 
ing strength at k ~ 1 h Mpc - 1 is probably the reason why 
iRudd, Zentner. & Krav tsov ( 2008j) did not detect it in their simula- 
tions. The box size that these authors used is relatively small, so that 
at those scales they were probably still plagued by finite volume ef- 
fects. Fo r the same reason , sim ulations using a l a rge bo x size, like 
those of Ijing et all d2006l) and Ivan Daalen et al.l d201 ll) did detect 
this feature, having the same magnitude as here. Despite its very 
small entity, the shape of this deviation in itself is interesting. If 
one looks at the effect of primordial non-Gaussianity on the matter 
power spect rum will find a very similar deviation if the skewnes s 
is negative jFedeli & Moscardinill2fJTrj| ; |Srnith & Markovidl201 lb . 
The scale of the deviation due to non-Gaussianity is maybe a little 
bit lower at z = 0, but most likely has a different redshift depen- 
dence than the scale of the one due to baryonic physics. A variation 
of ~ 1% as the one induced by gas physics could easily be mistaken 

with a primordial non-Gaussianity at the level of / NL 50, if one 

ignores the information at k > 2-3/7 Mpc -1 . Therefore, partic- 
ular care must be taken when interpreting matter power spectrum 
results in cosmological terms. 

Next, we turn to the quintessence cosmologies. In Figure [8] 
we show the same as in Figure [7j for each of the four dynamical 
DE models considered in the present work. The first thing that we 
note is that the qualitative shape of the total matter power spec- 
trum is the same as in the reference ACDM cosmology. On the 
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Figure 8. The total matter power spectra divided by the dark matter power spectra in the dark matter only simulations (cyan dot-dashed lines). For reference, 
we also show the dark matter only power spectra (black solid lines) and the total power spectrum for the fiducial cosmology, as in Figure[7](yellow solid line). 
Each panel refers to z = 0, and to a different quintessence cosmology, as labeled in the plots. Note that the numbers on the ordinate axis express a percentage. 



quantitative level, it looks like the depression at intermediate scales 
is slightly more pronounced and shifted at smaller scales for the 
SUGRA model as compared to the other quintessence models. This 
would make sense, because the SUGRA model is also the model 
where the gas cooling and star formation is more suppressed, hence 
it is expected to have a larger fraction of gas in the hot phase and 
stars concentrated at smaller scales. However, the differences are 
really minuscule, and they overall imply that the distribution and 
physics of the large scale hot gas are not significantly influenced by 
the quintessence model. 



7 DISCUSSION AND CONCLUSIONS 

We studied the correlation properties of dark matter, gas, and stars 
in cosmological models dominated by a dynamical DE component, 
making use of a suite of large cosmological simulations. The main 
results of our study can be summarized as follows. 

• In the fiducial ACDM cosmology the gas is much less clus- 
tered than the dark matter, more so for small scales, as expected for 



a part of it to cool down and condense into stars, as well as because 
of shock heating and gas pressure. The stars are in turn much more 
clustered than the other matter components, because they form in- 
side biased dark matter halos. The backreaction of gas onto dark 
matter is such that the clustering strength of the latter is increased 
by ~ 20% at small scales (k ~ 10 h Mpc -1 ). The total matter fol- 
lows a trend very similar (although not identical) to the dark matter. 

• The results of the previous point are in broad overall agree- 
ment with previous simulations of structure formation in pres- 
ence of various gas physical processes not including AGN feed- 
back. Nevertheless, there are a few quantitative differences. Among 
these, we find a decrement of the total matter power spectrum at in- 
termediate scales k ~ 1 h Mpc -1 , that some authors find and others 
do not. The reason for this is probably that the decrement is very 
small, ~ 1%, at the same order of our own modeling uncertainties. 
Also, the amount of the increment of clustering strength in the total 
matter power spectrum at small scales that we find is of ~ 20-30%, 
intermediate amongst other findings. 

• We used the physically motivated halo model in order to re- 
produce the backreaction of baryonic physics on the dark matter 
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power spectrum. In particular, we found that by increasing the av- 
erage concentration of dark matter halos by 17% irrespective of the 
mass accounts very well for the small-scale increase in dark mat- 
ter power due to the inclusion of gas and stars. By using this very 
simple prescription, the ratio of the dark matter power spectrum to 
its counterpart in the dark matter only simulation is reproduced at 
better than 1% at all scales of interest. 

• The ratios between various matter power spectra in 
quintessence cosmologies and their counterparts in the fiducial 
model have a very specific form that does not depend qualitatively 
on the precise cosmological model and the matter component. Par- 
ticularly, these ratios are flat at large scales, have a depression at 
k ~ 1 h Mpc~', and then raise again at small scales. At least part 
of this trend is due to the different erg values adopted in the various 
models in order to compensate the variations in the growth history. 

• By blindly applying the halo model to the quintessence cos- 
mologies, we manage to reproduce the overall shape of the dark 
matter power spectrum (normalized by the dark matter only spec- 
trum). However there are a few significant quantitative differences. 
One of them is that the mass function normalization in these mod- 
els ne eds to be artificia l ly low ered by a few percent with respect 
to the ISheth & Tormenl d2002h prediction, in order to match the 
dark matter power spectrum behavior at intermediate scales. This 
would mean that slightly less matter is locked into halos in the 
quintessence models, in agreement with a less efficient structure 
formation (see below). 

• Quantitatively, the power spectra of the various matter com- 
ponents hint at the star formation process being less efficient in 
quintessence models as compared to the ACDM cosmology, with 
the SUGRA model being the least efficient of all. For instance, the 
gas clustering strength in dynamical DE cosmologies is less de- 
creased, compared to the dark matter one, with respect to the fidu- 
cial model, implying that less gas condenses into stars. At the same 
time, the stars gain more power, implying that they can form only 
in more massive (and hence more biased) halos. The whole struc- 
ture formation process seem s to be less effective in quintessence 
cosmologies, as suggested bv lDe Boni et alj J20 10h . 

• As a further confirmation of this, we also found that the back- 
reaction of the gas onto the dark matter distribution is less severe 
for the quintessence cosmologies. As it turns out by using the halo 
model, the average concentration of dark matter halos needs to be 
increased by ~ 17% at z = in the ACDM model (see above), and 
only by ~ 10% in quintessence models. The exception to this is the 
RP model, which is the one with least dissimilarities with respect 
to the fiducial cosmology. 

• Finally, we also briefly studied the total matter power spec- 
trum, being the subject of cosmic shear observations. When com- 
pared to the dark matter only power spectrum, it displays a percent- 
level reduction of power at intermediate scales, driven by the gas 
component, and a sharp increase at small scales, driven by the stel- 
lar component. Contrary to the dark matter spectrum, the specific 
behavior of the total matter power cannot be accurately modeled 
with simple modifications to the halo model. 

As noted very recently bv lvan Daalen et al | djfjb.AGN feed- 
back, which was not included in our simulations, might have a sub- 
stantial impact on the large scale distribution of matter. Particu- 
larly, this has the effect of removing large amounts of gas from the 
centers of clusters and groups of galaxies, thus inhibiting the to- 
tal matter clustering at intermediate/small scales. The consequence 
is that the total matter power spectrum would be suppressed more 
and at larger scales compared to what we have shown in Figure [7] 



This also hints to the fact that differences with works where this de- 
pression is not noted might depend on the efficiency of the energy 
feedback from supernovae. The effect of AGN feedback would de- 
pend on the co-evolution of supermassive black holes with galax- 
ies, which would likely depend in turn on the expansion history 
of the Universe. It is therefore not straightforward to understand 
the effect of AGNs in dyna mical DE cosmo logies without further 
simulations. Nevertheless, van Daalen Tt al.1 d201ll) show that, for 
instance, by changing the normalization <x 8 by a significant amount 
has a very little effect on AGN feedback. 

While this paper was being finalized, a preprint with rel- 
evance to the p r esent work appeared on the archive, namely 
ISemboloni etail fcOllh . There the authors come up with 
a modification to the halo mod el similar to the one of 
Guillct, Tevssier. & Colombi d2010h , specifically overlapping to 
the dark matter density profile a /?- model for the gas density profile 
and a point mass representing the central star concentration. They 
show that in this way they manage to roughly reproduce the total 
matter power spectrum in units of the dark matter only one, with an 
accuracy of ~ 10%. Although their work is not yet complete, in that 
it lacks the backreaction of gas onto dark matter and the smoothly 
distributed gas component, it is an important step forward in mod- 
eling the large scale matter distribution on the Universe. 

Despite missing the AGN feedback ingredient, we believe the 
present work to be a substantial improvement in understanding the 
large scale gas and star distributions, as well as their backreaction 
on the dark matter power spectrum. All of these ingredient are go- 
ing to be of fundamental importance in order to construct a reli- 
able semi-analytic model describing the correlation function of all 
matter in the Universe. We also provide crucial information about 
the main differences in structure growth that arise when the back- 
ground expansion is dominated by quintessence. This is important 
in order to better understand the differences in the large scale mat- 
ter distribution that arise due to physically motivated quintessence 
scenarios, and hence to devise new tests for detecting a possible 
redshift evolution of the DE equation of state parameter. 
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